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Abstract. We study decoherence of two non-interacting qubits. The 
environment and its interaction with the qubits are modelled by random 
matrices. Decoherence, measured in terms of purity, is calculated in linear 
response approximation. Monte Carlo simulations illustrate the validity of this 
approximation and of its extension by exponentiation. The results up to this 
point are also used to study one qubit decoherence. Purity decay of entangled 
and product states are qualitatively similar though for the latter case it is slower. 
Numerical studies for a Bell pair as initial state reveal a one to one correspondence 
between its decoherence and its internal entanglement decay. For strong and 
intermediate coupling to the environment this correspondence agrees with the 
one for Werner states. In the limit of a large environment the evolution induces a 
unital channel in the two qubits, providing a partial explanation for the relation 
above. 



PACS numbers: 03.65.Yz,03.65.-w,03.65.Ud 

1. Introduction 

In recent experiments [Tl[2l[3j|4] it has been demonstrated that it is possible to protect 
ever larger entangled quantum systems, often arrays of qubits, ever more efhciently 
from decoherence. This improved protection of states makes a more detailed analysis 
of decoherence desirable. In the context of fidelity decay, a random matrix description 
[51 [B] is accessible and very effective in describing experiments [3 [HI [S]- A close 
connection between the dynamics of fidelity decay and decoherence has been shown 
in some instances [THl [H], which suggests to apply methods successful in one field 
to the other. Based on some previous work we shall therefore analyse in detail the 
consequences of using random matrix theory (RMT) to model decoherence [6l [121 US] 
for the case of two non-interacting qubits or a single qubit as the central system. 

Two perspectives make such a random matrix treatment particularly attractive. 
First, reduction of decoherence may, in some instances, be achieved by isolating some 
"far" environment (including spontaneous decay) to a degree that it can, to first 
approximation, be neglected. Then it can happen that the Heisenberg time of the 
relevant "near" environment is finite on the time scale of decoherence. In such a case 
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it becomes relevant that RMT shows, in hnear response approximation, a transition 
from hnear to quadratic decay at times of the order of the Heisenberg time. This 
behaviour is seen with spin chain environments |14| , and is essential for the success of 
the theory in describing the above mentioned experiments of fidelity decay. Note also 
that the concept of a two stage environment has been used for basic considerations 
[15] . Second, the long term goal must be to describe in one theory the decay of fidelity 
that includes undesirable deviations of the internal Hamiltonian of the central system, 
together with decoherence. 

This study is focussed on weak coupling of one or two qubits to an environment, 
and thus we can use the correlation function approach proposed for purity decay in 
echo-dynamics |16j , treating the coupling as the perturbation [6] . The linear response 
approximation will be sufficient and in this approximation the ensemble averages, 
which we have to take in any RMT model, are feasible though somewhat tedious. 
Exact solutions, which exist in some instances for the decay of the fidelity amplitude 
[T71[TH], seem to be out of reach at present, because they would require the evaluation 
of four-point functions. 

Assume that the qubit(s) are initially in a pure state, and evolve under their 
own local Hamiltonians. The qubit(s) are coupled independently via random matrices 
to a large environment in turn described by another random matrix. The coupling 
to the environment gives rise to decoherence. Averaging both the coupling and 
the environment Hamiltonian over the RMT ensembles yields the generic behaviour. 
We then focus on the correspondence between decoherence and entanglement decay 
for a Bell state. We find a relation between this two quantities, that can be 
partially understood assuming unitality. This relation in most cases coincides with the 
corresponding one obtained for Werner states. For that part, it is essential that the 
two qubits do not interact. Otherwise, the coupling between the qubits would act as an 
additional sink (or source) for internal entanglement - a complication we wish to avoid. 
However, a residual coupling between the qubits should be taken into account at some 
point [19]. A preliminary study of some aspects of this system has been presented 
in [13], but here we shall present a more general picture: First, we allow for local one- 
qubit Hamiltonians. They have a considerable effect on decoherence. Second, we cover 
exhaustively aZZ possible initial states. Third, we study the consequences of the whole 
system being time- reversal invariant. Indeed, we present the mathematical procedure 
in detail using both the Gaussian unitary (GUE) and the Gaussian orthogonal (GOE) 
ensembles [2Qj [21] for the description of the environment and the coupling. The 
two ensembles correspond to time reversal invariance (TRI) breaking and conserving 
dynamics respectively. All these cases are treated on equal footing. 

In section [2] we shall state the model, recall the linear response formalism for 
echo dynamics, and show how it can be adapted to forward evolution. In section [3] 
we shall discuss decoherence of a single qubit, which already shows some interesting 
features by itself and serves as a warm up for the more complicated case of two 
qubits. The details of the calculations are given in the appendices. This shall 
be treated in section IH starting from the situation where one qubit is a spectator, 
exerting influence only through entanglement. More general cases are analysed and 
reduced to the spectator situation. Most analytic results will be accompanied by 
Monte Carlo simulations. In section [5] we investigate the correspondence between 
decoherence (purity) and entanglement (concurrence) as both evolve in time. We end 
with conclusions in section [6] 
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2. General considerations 

In this section, we describe the general framework for our study of decoherence. We 
consider one- and two-qubit systems, coupled to an environment which is modelled 
using random matrix theory. The one-qubit case is introduced in section [2.11 the two- 
qubit case in section [2T2l The quantities studied in this paper are defined in scction [2T3l 
and the tools used to solve the problem are introduced in section [2^ 

2.1. The model for one qubit 

We describe decoherence by considering explicitly the additional degrees of freedom 
(henceforth called "environment") which are interacting with the qubit. We follow 
the unitary time evolution of a pure, initially separable, state in the product space 
T-L = Til ® Tie, where Tii (of dimension two) and TL^ (of dimension N^) denote the 
Hilbert spaces of the qubit and the environment, respectively. As time evolves, the 
qubit and the environment become more and more entangled, which means that after 
tracing out the environmental degrees of freedom, the state of the qubit becomes more 
and more mixed. For such a setup, the Hamiltonian is of the following form 

Hx=Hi®U + li®H^ + AT/i,c = Hi+H^ + XVi,^ . (1) 

Here, Hi represents the Hamiltonian acting on the qubit, the Hamiltonian of the 
environment, and T^.e the coupling between the qubit and the environment. The real 
parameter A controls the strength of the coupling. For notational ease, and if there is 
no danger of confusion, we omit extensions with the identity, as it is done in ([1]). 

We describe both the coupling and the dynamics in the environment within 
random matrix theory. To this end. He and Vc^i are chosen from either the GUE or the 
GOE, depending on whether we wish to describe a TRI breaking or TRI conserving 
situation. The Hamiltonian Hi implies another free parameter of the model, namely 
the level splitting A of the two level system representing the qubit. 

We shall study the time evolution of an initially pure and separable state 

m = o)> = ivi) ® iVe) , (2) 

where jV'i) € H-i and \il}c) S He. At any time t, the state of the whole system is thus 
\ip{t)) = exp(— iiiJA)IV'(O)), and the state of the single qubit is trc |-0(i))(V'(OI- While 
the state of the qubit \tpi) is yet another free parameter in our model, we assume 
the state of the environment j-i/jc) to be random. This means that the state is chosen 
from an ensemble which is invariant under unitary transformations. In practice, this 
means that the coefficients are chosen as complex random Gaussian variables, and 
subsequently the state is normalized. 

2.2. The model for two qubits 

For the two qubit case, we need to introduce the additional Hilbert spaces H2 of 
dimension two and He' of dimension Nc' , which correspond to the second qubit and 
its environment respectively. As in the one-qubit case, the dynamics in the space of 
the two qubits Hi®H2 (henceforth called "central system") is obtained by considering 
the unitary evolution in whole Hilbert space Hi®H2®Hc® He' and subsequently 
tracing-out the environmental degrees of freedom (henceforth called "environment" ) . 
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As initial states, we use pure states in the whole Hilbert space, which are separable 
with respect to any combination of Hi ® Ti.2, He, and He'- 

1^(0)) = |V'l2)|^e)|^e'), |V'12> eHl<^H2, IV'c) GHe, andl^e') eT^e' . (3) 

Typically, we assume that at t = the two qubits (and hence the state \1p12)) to be 
entangled. This is the main new ingredient, as compared to the one-qubit case. In 
this work we shall consider three different dynamical scenarios, all explicitly excluding 
any interaction between the two qubits: 




Figure 1. Schematic representations of the different dynamical configurations 
studied in this article: (a) the spectator Hamiltonian, (b) the separate, and (c) 
the joint environment Hamiltonian. 



(a) The spectator Hamiltonian: In this case the two qubits evolve under separate 
Hamiltonians Hi and H2 acting on the Hilbert spaces Tii and 7^2, respectively. 
In addition, we assume that only the first qubit is coupled to an environment, 
such that the total Hamiltonian reads 

Hx = Hi+H2 + H, + AFi,e , (4) 

where A, Hi, H^ and Vi^o are defined as in ([1]). This situation is shown 
schematically in figure HJa) . If we choose an initial state where the two qubits 
are already entangled, this provides the simplest situation which allows to study 
entanglement decay. If the two qubits are initially not entangled, the process 
reduces effectively to the one-qubit case, described previously section 12.11 The 
special case Hi = H2 = has been considered in Ref. [T5] . 

(b) The separate environment Hamiltonian: This case differs from the spectator 
model only inasmuch as the second qubit is also coupled to an environment. 
Both environments are assumed to be non-interacting. The total Hamiltonian 
then reads 

Hx,M = Hi+H2+H, + ffe' + Ai T/i,e + X2 T^2,c' , (5) 

where V2.0' and H^i describe the coupling to - and the dynamics in the additional 
environment. Both quantities are chosen independently from the respective 
random matrix ensembles, in perfect analogy with Vi^c and H^. The real 
parameters Ai and A2 fix the coupling strengths to either environment. This 
model, see figure [Ub), may describe two qubits that are ready to perform a 
distant teleportation, where each of them is interacting only with its immediate 
surroundings. It can also represent a pair of qubits that, although close together, 
interact with different and independent degrees of freedom. 
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(c) The joint environment Hamiltonian: The third case, shown in figure [He), 
describes a situation in which both qubits are coupled to the same environment, 
even though the coupHng matrices are still independent. In that case, the total 
Hamiltonian reads 

Hx,,x, =H,+H2 + H, + X, Vi^, + X2 , (6) 

where V2,o describes the coupling of the second qubit to the environment. It is 
chosen independently from the same random matrix ensemble as Vi^o- 

2.3. The measures 

As a measure of decoherence we use purity. This measures the degree of mixedness of 
a density matrix p representing a physical system. Purity is defined as 

P[p\^trp\ (7) 

It reaches a maximum value of one when the state is pure, i.e. when p = \ip){4'\ for 
some I '!/')■ Otherwise it is less than one and reaches a minimum when the state is 
completely mixed. We use purity instead of e.g. the von Neumann entropy because 
it is simpler to handle from an algebraic point of view. Note that for one qubit both 
quantities are equivalent, and for two qubits they contain similar information [22j . 

The other measure considered here is concurrence. It quantifies the degree of 
entanglement of two qubits in a pure or mixed state. It is a measure used extensively 
in the literature [53] that is straightforward to calculate. It is closely related to the 
entanglement of formation, which measures the minimum number of Bell pairs needed 
to create an ensemble of pure states representing the state to be studied [24]. Given 
a density matrix p representing the state of two qubits, concurrence is defined as 

C{p) = max{0, Ai - A2 - A3 - A4} (8) 

where A^ are the eigenvalues of the matrix \J pipy ® <^y)p*{<^y ® fj,) in non-increasing 
order. The superscript * denotes complex conjugation in the computational basis and 
CTj, is a Pauli matrix [25j . The concurrence has a maximum value of one for Bell states, 
and a minimum value of zero for separable states. Furthermore, it is invariant under 
bilocal unitary operations. 

Both, purity and concurrence, provide a measure of entanglement, but in the 
present context they quantify completely different aspects of the problem. Purity will 
be used to quantify the entanglement between the pair of qubits and the environment 
i.e. decoherence. Concurrence will be used to measure the entanglement within the 
pair. Although they are not fully independent, information about one provides little 
knowledge about the other. E.g. a state with purity one can have full entanglement 
(a Bell pair) or null entanglement (a separable pure state). 

2.4. Echo dynamics and linear response theory 

We shall calculate the value of purity as a function of time analytically, in a 
perturbative approximation. To perform this task it is useful to consider each of 
the above Hamiltonians [([J), ([4]), (O, and ([6])], as composed by an unperturbed part 
Ho and a perturbation XV with A a parameter. The unperturbed part corresponds 
to the operators that act on each individual subspace alone whereas the perturbation 
corresponds to the coupling among the different subspaces; for example in the one 
qubit case Hq = He + Hi and V = 14,1- For the case of two qubits, A can be taken 
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as maxjAi, A2}. We shall use the tools developed in [Appendix A| for a linear response 
formalism in echo dynamics. In order to do so we must formally write the problem in 
this language. 

We write the Hamiltonian as 

H^=Ho + XV, (9) 

and introduce the evolution operator and the echo operator defined by 

Ux{t) = e-'^-* , Ahit) = Uoit)^Uxit) , (10) 

respectively {h = 1). For the calculation of purity at a given time t, we replace 
the forward evolution operator U\ hy the corresponding echo operator M\. Even 
though the resulting states are different, i.e. p{t) = tr^.c' U\{t)pU\{t) ^ p^'^ {t) = 
ivc.c' M\{t)pM\{t)^ they are still related by the local (in the two qubits and the 
environment) unitary transformation C/o(i). Since such transformations do not change 
the entanglement properties, it holds that P{t) — P[p{t)] = P[p'^^{t)]. This step is 
crucial, since the echo operator admits a series expansion with much larger range of 
validity (both, in time and perturbation strength). The numerical simulations are all 
done with forward evolution alone as they require less computational effort. 
The Born expansion of the echo operator up to second order reads 

Mx{t) = 1 - iXI{t) - X^J{t) + 0(A3) , (11) 

with 

I{t) = f drVir), J{t) ^ [ At I dT'V{T)V{T') (12) 
Jo Jo Ja 

and V{t) = Uo{tyVUoit) being the coupling in the interaction picture. Using this 
expansion we calculate the purity of the central system, averaged over the coupling 
and the Hamiltonian of the environment. 



3. One qubit decoherence 

We first study the GUE case with and without an internal Hamiltonian governing the 
qubit. The next step is to work out the GOE case. There we concentrate on the 
case with no internal Hamiltonian governing in the qubit since we want to keep the 
discussion as simple as possible to focus on the consequences of the weaker invariance 
properties of the ensemble. 

Here, we study the evolution under the Hamiltonian 

Hx = Hi+H, + X yi,e = Ho + X Fi,e , (13) 

where the initial state 

^'o = |^l)(V'l|®IV'c)(V'c| , (14) 

is the product of a fixed pure state jV'i) G TLi and a random pure state \i>c) G Ti-c'i see 
([2]). At any later time, the state of the qubit and its purity are given by 

p{t) = tre Ux{t) go Ux{t)^ P{t) = tri p{t)^ . (15) 

In [Appendix A[ we compute the average purity {P{t)) as a function of time in the 
linear response approximation (jlip . following the steps outlined in section [2^ The 
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average is taken with respect to the coupUng Vi.c, the random initial state jV'e) s^^id 
the spectrum ol H^. In the limit of oo, we obtain [ (|A.7|) . (|A.10|) ] 

{P{t)) = l-2\^ J dT J dr'ReAji(r,r') + 0(A*) , (16) 

A.n{T, t') = [Ci(|r - r'|) - S,{t - t')]C{\t - r'|) + xgoe[1 - S[{-t ~ r')] , (17) 

where xgoe = 1 for the TRI case, and xgoe = for the non-TRl case. The correlation 
functions Ci(t), S'i(r), S'J(t), and C{t) are defined in [Appendix B[ C{t) deserves 
special attention, since all the dependence of the environment is via the function 



1 



Ef=ie-'^^* )=C{t) = l + S{t/TH)-brit/TH), (18) 



where the Ej's are the eigenenergies of He and th is the corresponding Heisenberg 
time. The two-point form factor, b2^\T), is known analytically for the GUE (/3 = 2) 
and the GOE (/? = 1), which are the two cases treated here ^21) . 

3.1. The GUE case 

We are now in the position to give an explicit formula for (P(t)) in the GUE case. 
This formula will generally depend on some properties of the initial condition I'f/'i). 
We wish to write lipi) in the most general way, but grouping together cases, equivalent 
due to the invariance properties of the problem. 

Recall that H = + Hi + XV represents an ensemble of Hamiltonians in which 
He G GUE and V G GUE, whereas Hi (together with the initial condition \ipi)) 
remains fixed throughout the calculation. The operations under which the ensemble 
is invariant are local (with respect to the partitioning of the Hilbert space into Hi and 
He), unitary (due to the invariance properties of the GUE), and leave Hi invariant. 
Hence the transformation matrices must be of the form U ® exp(ia-ffi) with a a real 
number and U a unitary operator acting on He- 

This freedom allows to choose a convenient basis to solve the problem. On one 
hand, it allows to write H^ in diagonal form (as done in [Appendix Ap , and on the 
other hand, we can use it to represent the initial state of the qubit in such a way that 
there is no phase shift between the two components of the qubit. This can be achieved 
by appropriately choosing a (see also the discussion in section |4T|) . We thus write, 
without loosing generality 

=cos0|O) +sin0|l) , (19) 

where |0) and |1) denote the eigenstates of Hi. Notice that if </) = and </> = tt/2, {ipi 
is an eigenstate of Hi . Finally, we choose the origin of the energy scale in such a way 
that the Hamiltonian of the qubit can be written as Hi = (A/2)|0)(0| - (A/2)|l)(l|. 
Hence, is the only relevant parameter describing the initial state. 

We obtain the average purity from the general expression in (|16p and ^Tl} . For a 
pure initial state pi = \'ipi){ipi\ the relevant correlation functions ReCi(T), Si{t) and 
C{t) are given in (|B.10|1 . (jB.12p . and (jB.Sp . respectively. Using the symmetry of the 
resulting integrand with respect to the exchange of t and r', we find 

(P(i))-1-4A2 f dT Tdr' C(t') [1-50 (1-cosAt')] +0(A4,iV-i) (20) 
Jo Jo 



with 



A A 3 + cos(4(A) 
cos^ (j) + sin"* (j) = (21) 
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quantifying the "distance" between and the eigenbasis of Hi . 

Let us consider following two limits for Hi. The "degenerate limit", where the 
level splitting A is much smaller than the mean level spacing do = 2'k/th of the 
environmental Hamiltonian, and the "fast limit" , where the level splitting is much 
larger. In the latter case, the internal evolution of the qubit is fast compared with the 
evolution in the environment. 

The degenerate limit leads to the known formula (13j 

FoW^l-AVrnW, (22) 



with 



2t^ 

2tTii + — ifO<t<TH, 
3th 



2t'^ + if i > th 

3 

The result does not depend on the initial state of the qubit. Due to the degeneracy all 
states are eigenstates of Hi and thus equivalent. The leading term of the purity decay 
is linear before the Heisenberg time and quadratic after the Heisenberg time. Similar 
features were already observed in fidelity decay and purity decay in other contexts [5] . 

In the fast limit (A ^ do), purity is obtained from ^U\\ by replacing cos Ar' with 
one when it is multiplied with the S function [see (jlSp j. and with zero everywhere else. 
For finite A'e care must be taken, since we are assuming Zeno time (which is given 
by the "width of the (5- function" ) to be much smaller than all other time scales, such 
that A <C iVo dc- The resulting expression is 

PAt) - 1 - A2[(1 - g^)frAt) + 2g^tm] . (24) 

Typically (depending on the initial state), this formula again displays a dominantly 
linear decay below the Heisenberg time, and a dominantly quadratic decay above, 
similar to (HH). 

In figure [2] we compare numerical simulations of the average purity {P{t)) 
(symbols) with the corresponding linear response result (dashed lines) based on ([^0]) . 
The numerical results are obtained from Monte Carlo simulations with 15 different 
Hamiltonians and 15 different initial conditions for each Hamiltonians. We wish to 
underline two aspects. First, the energy splitting in general leads to an attenuation 
of purity decay. Even though a strict inequality only holds for the limiting cases, 
Ppit) > Pr>{t) (for t 7^ 0), we may still say that increasing A tends to slow down 
purity decay. This result is in agreement with earlier findings on the stability of 
quantum dynamics |26| . Second, for the fast limit and an eigenstate of Hi [g^ — 1) 
we find linear decay even beyond the Heisenberg time. A similar behaviour has been 
obtained in [B], but there an eigenstate of the whole Hamiltonian was required. 

In [13] it was shown that exponentiation of the linear response result leads to 
very good agreement beyond the validity of the original approximation. We use the 
formula ((CT|) 

1-Plr(<)" 



PELR(t) = Poo + (l-Poo)exp 



l-Po< 



(25) 



where Pi,B.{t) is truncation to second order in A of the expansion (|20|) . and Poo = 1/2 
the estimated asymptotic value of purity for t — s- oo, see |Appendix C[ From figureOwe 
see that the exponentiation indeed increases the accuracy of the bare linear response 
approximation. 




Figure 2. Numerical simulations for the average purity of one qubit using non- 
TRI Hamiltonians, as a function of time in units of the Heisenberg time th of the 
environment. For the coupling strength A = 0.01 and A^c = 2048, we show the 
dependence of (P(t)) on A (the level splitting) and (j> (characterizing the initial 
state) of the internal qubit: A = 0, (/i = 7r/4 (green circles), A = 8,<^ = 7r/4 (blue 
squares), and A = 8,0 = (red triangles). The corresponding linear response 
results (dashed lines) and exponentiated linear response results (solid lines) are 
based on 1201 1 and I I25I I. where Poo is given in [Appendix C| Note that the level 
splitting A tends to slow down decoherence. 



3.2. The GOE case 

We now drop Hi leaving Hq = He, resulting in 

Hx = H, + XV (26) 

where H^ £ GOE acts on He and V e GOE on Hc^^Hi. The resulting ensemble of 
Hamiltonians is invariant under local orthogonal transformations. In the environment, 
this allows again to diagonalize H^. In the qubit it allows rotations of the kind 
exp(iatTy) g 0(2). If such transformations are represented on the Bloch sphere, they 
become rotations around the y axis. Hence, they can take any point on the Bloch 
sphere onto the ccy-plane. Supposing this point represents the initial state, it shows 
that we may assume jV'i) to be of the form 

I*,) ^ Ml , ,„) 

In this expression, 7 e [— 7r/2,7r/2] denotes the angle of the point representing the 
initial state with the xz-plane (see figure [3]). 

In order to obtain the linear response expression for {P{t)) we again make use of 
(|16p and (|17p . However, apart from the correlation functions used in the GUE case, 
we have now to consider in addition S[ (r) , as given in (jB.17p . The special case Hi ~ 
can simply be obtained by setting A = 0. This yields 

^ji(r,r')-C(|T-r'|) + sinS- (28) 
After evaluating the double integral in (fT6|) . we obtain 

(P(t)) = 1 - A2 {t2 [3 - cos(27)] + 2tTH - 25^^^ (t) J , (29) 
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Figure 3. Any pure initial state of the qubit can be mapped onto the Bloch 
sphere. Here, we show the angle 7 defined in ||27|I in colour code. Regions of 
a given colour represent subspaces which arc invariant under the transformation 
exp(icKTy). 

where 

Bi'\t)^2 fdr r dr' bi'\T' /th) (30) 
Jo Jo 

is the double integral of the form factor. It can be computed analytically, but the 
resulting expression is very involved |5j. For our purpose it is sufficient to note that 
for t <C th, B2^\t) cx t^ (as in the GUE case), whereas for t^ th, t — B2^\t) grows 
only logarithmically. 




1 2 3 4 5 6 

t/TH 



Figure 4. The averages of the purity and the von Neumann entropy for initial 
states, "located" in different regions on the Bloch sphere. These regions are 
characterized by the angle 7 = (green squares) and 7 = 7r/2 (blue circles). The 
lines corresponding to the formula (I29I I with the appropriate value of 7. Here, 
TVe = 2048 and A = 10"^. 

In figure[4]we show {P{t)) for 7 = (green squares) and for 7 ~ tt/2 (blue circles). 
In contrast to the GUE case (degenerate limit), the average purity depends on the 
initial state via the angle 7. The fastest decay of purity is observed for 7 = 7r/2, where 
the image under the time reversal operation becomes orthogonal to the initial state. 
The slowest decay is observed for 7 = 0, which characterizes states which remain 
unchanged under the time reversal symmetry operation. These statements can be 
directly translated to the von Neumann entropy S. For one qubit it has a one to one 
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relation with purity 

S{P) = h ( ii^^) + h ( ^^4^ ) > M-) = -xlofo ^. (31) 



Observe the entropy scale on figure ID 

Due to the dependence of (|29p on 7, different initial conditions in the Bloch sphere 
will yield different behaviours of purity. However for a fixed value of 7 we numerically 
check that there is self averaging (not shown) , meaning that different members of the 
ensemble behave the same as the average, in the large limit. As for one qubit the 
von Neumann entropy is a function of purity, this behaviour will also be observed in 
this quantity. The consequences of the weaker invariance properties of the GOE will 
be analysed in detail in a more general framework in a latter paper. 

4. Two qubit decoherence 

In this section, we address the question whether entanglement within a given system 
affects its decoherence, if the system is coupled to an environment. We discuss the 
three different situations described in section [ 



4.I. The spectator Hamiltonian 

The spectator Hamiltonian discussed in section [^21 reads 

Hx = Hi+ H2 + H^, + AFc,i . (32) 
In order to study the decoherence of the initial state 

|V'l2)(^12|<^|^c)(V'c| , (33) 

we use the fact that the dynamics in 7^2 also decouples from that in Tii and Tic, 
that H2 commutes with all other terms in the Hamiltonian. The quantum echo of go 
after time t is 

g^' (t) - ( I2 ® Ma (t)) go{l2® Ml (t) ) . (34) 
Since g^ (t) remains a pure state in Hi 7^2 (X" He for all times, 

F(t) = trpc(t)' =trpe(t)' pcit) = tr, g^^ (t) p,{t) = tr, g^\t) . (35) 
As the echo operator acts as the identity on the second qubit, 

p,{t) = tn Mx{t){tv2 go)Ml{t) = tn Mx{t) {pi (g, \^c){^c\ ) Ml{t) , (36) 

where pi — tr2 |V'i2)(V'i2|- We may therefore compute the purity of the spectator 
model, without ever referring explicitly to the second qubit. Any dependence of the 
decay of the purity on the entanglement between the two qubits is encoded into the 
initial density matrix pi. This also implies that we can use the results obtained 
in [Appendix A[ and hence ITBl) and ITTl) remain valid. The only difference is that for 
the correlation functions Ci(t), 6*1 (t), and S[{t), we now have to insert the respective 
expressions which apply for mixed initial states of the first qubit. These expressions 
are given in Appendix B[ 
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|0) = IO2 




Figure 5. We see a plot to help visualize tlie way the initial condition 
is parametrized. On the left, qubit one, has an internal Hamiltonian. Its 
eigenvectors (|0) and |1)) are represented in blue. The z axis is chosen parallel to 
the vector |0). The x axis is chosen in such a way as to make both |0i) and |li) 
have real coefficients i.e. such that the xz plane contains |0i) and |li). On the 
right we represent the second qubit where we have absolute freedom to choose the 
basis (even if an internal Hamiltonian is present) , and thus we choose it according 
to the natural Schmidt decomposition. 



4.I.I. The GUE case: We again wish to write the initial condition in its simplest 
form. We must respect the structure of ([U, but take advantage of all its invariance 
properties. Given a fixed i/i, the ensemble of Hamiltonians is invariant under local 
operations of the form [/Ar_,(K)expiai/i(g)f72 where Um^ S U{Nc) is any unitary operator 
acting on the environment, a a real number, and U2 G ^(2) is any unitary operator 
acting on the second qubit. 

The freedom in the environment allows to write in diagonal from, whereas the 
freedom within the qubits allows to choose a basis {|0), |1)} ® {|0), |1)} in which the 
state can be written as 

1^12) =cos6'(cos(^|0) +sin(/)|l))|0) +sin6'(sin0|O) -cos<?!)|l))|l),(37) 

and still. Hi — ■y|0)(0| — ■y|l)(l| is diagonal. To find this basis we start using the 
Schmidt decomposition to write 

|Vi2> =cos0|OiO2) +sin0|iii2) (38) 

with {|0i), |li)} being an orthonormal basis of particle i. For the first qubit, we fix the 
z axis of the Bloch sphere (containing both |0) and |1)) parallel to the eigenvectors 
of Hi, and the y axis perpendicular (in the Bloch sphere) to both the z axis and 
|0i). The states contained in the xz plane are then real superpositions of |0) and 
|1), which impHes that |0i) = cos0|O) + sin(/()|l) and |li) = sin(/)|0) — cos(^|l) for 
some (j). In the second qubit it is enough to set |0) = IO2) and |1) — II2). This 
freedom is also related to the fact that purity only depends on tr2 |'0i2)(^i2|- A 
visualization of this procedure is found in figure [51 The angle 6 £ [0, 7r/4] measures 
the entanglement (C(|'0i2)(V'i2|) — sin 26*) whereas the angle (p € [0, 7r/2] is related to 
an initial magnetization. 

The general solution for purity using this parametrization is 

P{t)^l-A\^ fdr r At'C{t') W^l+gfl cos Ar' ] + 0(A^ TV"!) , (39) 
Jo JO 
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where the geometric factors 5^^^ £ [0, 1/2], and g^^^ G [1/2, 1] are expressed as 

5e!i = 3e(l-30)+.90(l-3e) (40) 

g^'l ^ 2(1 ~gg)~g^{l~2gg), (41) 

in terms of the functions g^ and gg, defined in (|2ip . Both geometric factors are shown 
in figure [51 Equation ([59)) is obtained from (|16p and p7)) by insertion of the ([B.lOp 
and (|B.16P for ReCi(T) and 6*1 (r), respectively. 




Figure 6. Visualization of the geometric factors gg, g^^^, and g^^^ from left to 

right respectively. For g^^^ we see that for pure eigenstates of Hi its value is zero. 
This leads to a higher qualitative stability of this kind of states. 



We consider again two limits for A. In the degenerate limit (A <C I/th) purity 
decay is given by 

P^{t) = l-X^{2~gg)frAt), (42) 

where fmit) is defined in ([^51) . The result is independent of since a degenerate 
Hamiltonian is, in this context, equivalent to no Hamiltonian at all. The 0-dependence 
in this formula shows that an entangled qubit pair is more susceptible to decoherence 
than a separable one. 

In the fast limit (A 3> 1/th) we get 

(i) = 1 - [sg/xH {t) + 2rH5Si] • (43) 

For initial states chosen as eigenstates of Hi we find linear decay of purity both below 
and above Heisenberg time. In order for pi to be an eigenstate of Hi it must, first 
of all, be a pure state (in Hi). Therefore this behaviour can only occur ii 9 — or 
— tt/2. Apart from that particular case, we observe in both limits, the fast as well 
as the degenerate limit, the characteristic linear/quadratic behaviour before/after the 
Heisenberg time similar to the one qubit case. 

In figure [71 we show numerical simulations for {P{t)). We average over 30 
different Hamiltonians each probed with 45 different initial conditions. We contrast 
Bell states (0 = 7r/4, 9 = 7r/4) with separable states {(j) = 7r/4, 9 = 0), and also 
systems with a large level splitting (A — 8) in the first qubit with systems having 
a degenerate Hamiltonian (A = 0). The results presented in this figure show that 
entanglement generally enhances decoherence. This can be anticipated from figure [HI 
since for fixed 0, increasing the value of 9 (and hence entanglement) increases both 
gg^\ and gg'^l- At the same time wc find again that increasing A tends to reduce 




0.5 1 ij^^ 1.5 2 2.5 

Figure 7. Numerical simulations for the average purity as a function of time in 
units of the Hcisenberg time of the environment (spectator configuration, GUE 
case). For the coupling strength A = 0.03 we show the dependence of {P{t)) on 
the level splitting A in Hi and on the initial degree of entanglement between 
the two qubits (in all cases ip = n/ A): 8 = (separable states), A = 8 (black 
triangles); 6 = tt/4 (Bell states), A = 8 (blue rhombus); (9 = 0, A = (green 
squares); 6 = tt/4, A = (red circles). The corresponding linear response results 
(dashed lines and exponentiated linear response results (solid lines) are based on 
l|39| l and l|25|l . They are plotted with the same colour, as the respective numerical 
data. In all cases A^e = 1024. 

the rate of decoherence, while a strict inequaUty only holds among the two limiting 
cases (just as in the one qubit case). From g^^^ = 2 — g^^^ — gg, it follows that 

(Pp — Pd)/A^ = ffg^^i/rnC^) ^ 2iTH] > 0. Therefore, for fixed initial conditions and t 
greater than 0, Ppit) > Poit). This is the second aspect illustrated in figure [71 

In order to extend the formulae to longer times/smaller purities we exponentiate 
them using the results of |Appendix C[ The numerical simulations (see figure [7]) agree 
very well with that heuristic exponentiated linear response formula. In one case (blue 
rhombus) where the agreement is not so perfect, we found that it is the inaccurate 
estimate Poo = 1/4 of the asymptotic value of purity, which leads to the deviations. 

4-1.2. The GOE case: Let us consider the GOE average of both He and Vi^c- We 
have to take into consideration the weaker symmetry properties of the GOE ensemble 
as compared with the GUE. When averaging and Vi,c over the GOE, we are again 
confronted with the fact that the invariance group is considerably smaller than in 
the GUE case. In this context the initial entanglement between the two qubits has a 
crucial importance since it "transports" the invariance properties from the spectator 
to the coupled qubit. 

For the sake of simplicity, we focus on the degenerate limit setting Hi — Q. Note 
that on the basis of the results in |Appendix A[ the general case can be treated similarly 
and the corresponding result will be presented at the end of this subsection. 

We first specify the operations under which the spectator Hamiltonian ((32)) . 
considered as a random matrix ensemble, is invariant. As both, the internal 
Hamiltonian of the environment and the coupling, are selected from the GOE the 



Decoherence of two qubit systems: A random matrix description 



15 



invariance operations form the group 

0{Nc) X 0{2) xU{2), (44) 

and have the structure On^ <X) expijaay) (8) U2, with On^ being an orthogonal matrix 
(acting in Tic), ct a real number, and U2 a unitary operator acting on the spectator 
qubit. 

The direct product structure of the invariance group oblige us to respect the 
identity of each particle, but allows to analyse each qubit separately. For instance, if we 
would replace the random coupling matrix Vi_o with one, which involves both qubits, 
the invariance group would be 0{Nc) x 0(4). As a consequence, purity decay would 
become independent of the entanglement within the qubit pair: for any entangled 
state one can find a orthogonal matrix which maps the state onto a separable one. 

We write the initial condition \tjj12) as in For the coupled particle follows the 
same analysis made in section [3?2] We can thus write |0i) = 2~^/^[|0)+exp(i7)|l)] and, 
in order to respect orthogonality, |ii) = 2~^/^ exp(iC)[|0) — exp(i7)|l)]. For the second 
qubit we have the same complete freedom as in (|4.1.ip . We thus select IO2) = |0) and 
II2) = exp(— iC)|l) to erase the relative phase in the first qubit and finally write the 
initial state as 

cosg(|0)+e'^|l))|0)+sin^(|0)-e'^|l))|l) 

m2) = • (45) 

The average purity is still given by the double integral expression in (fT6|) . However, in 
the present case the mixed initial state pi = tr2 |V'i2)(V'i2| must be used. For A = 0, 
the resulting integrand reads 

Aji(r, r') = (2 - ge)Ci\T ~ r'|) + 1 - .ge + i^ge - 1) sin^ 7, (46) 

where Ci(t), S'i(r), and S[{t) are given in (jB.lOp . (|B.16p . and (|B.20p . respectively. 
Evaluating the double integral, we obtain 

(P(t)) = 1 - |t2[4_ 2cos2(26l)cos2 7] + (4- 2^9) trn - B!^^\t) | , (47) 

where B!^^\t) is given in (|30p . As in the GUE case, this result depends on the 
entanglement of the initial state, and, as in the one-qubit GOE case, it also depends 
on 7. Again it turns out that Bell states are more susceptible to decoherence than 
separable states. Note however, that purity as a function of 9 is not monotonous. 
Hence, a finite increase of entanglement does not guarantee that the purity decreases 
everywhere in time. For separable states, ge — 1, 'we retrieve formula (j29p . However, 
for completely entangled states, 9 = 7r/4, and the dependence on 7 is lost. This is 
understood from a physical point of view, noticing that any bilocal unital channel (in 
particular unitary operations) on a Bell state can be reduced to a single local unitary 
operation acting on a single qubit i.e. for any local unitary operation C/1,2, there exists 
U[ such that 

Ui (g) C/2|Bell) = U[(gi l|Bell) (48) 

(|Bell) is any 2 qubit pure state with C — 1, e.g. |00) + |11)). We can then say that 
the invariance properties in the second qubit are inherited from the first qubit via 
entanglement. 

Let us obtain the standard deviation for the different possible initial conditions in 
the qubits. We want to analyse the situation separately for a fixed value of concurrence. 
Then, as the invariant measure of the ensemble of initial conditions, and fixing the 



Decoherence of two qubit systems: A random matrix description 



16 



amount of entanglement, we shall use the tensor product of the invariant measures 
in each of the qubits. Since there is no dependence of (|T7)) on the second qubit, the 
appropriate invariant measure is trivially inherited from the invariant measure for a 
single qubit. The resulting value for the standard deviation is 

Based on [Appendix A| and [Appendix B| we can also obtain the average purity 
for A ^ 0. The parametrization of the initial states is more complicated since two 
preferred directions arise, one from the eigenvectors of the internal Hamiltonian and 
the other from the invariance group. The result can be expressed in the form given in 
(flB. with 



ReAn{r,T')^C{\T~r'\) U]l + gi'l cosA(r-r')] 

+5^ - (1 - & cos[A(t + r') - 2rj] + 0(A^ N'^) . (50) 

The angle t] is related to a phase shift between the components of any of the 
eigenvectors of the initial density matrix pi . 

4-2. The separate environment Hamiltonian 

We proceed to study purity decay with other configurations of the environment. 
Consider the separate environment configuration, pictured in figure [ijb). The 
corresponding uncoupled Hamiltonian is 

Ho = Hi+H2+H, + H,, (51) 

and the coupling is 

W = \iV,,i + \2V,,,2- (52) 

From now on we assume that the internal Hamiltonians of the environment and the 
couplings are chosen from the GUE. Generalization for the GOE can be obtained along 
the same lines using the corresponding results of section 14.1.21 The initial condition 
has a separable structure with respect to both environments, see ([3]). Next we calculate 
the coupling in the interaction picture. It separates into two parts acting on different 
subspaces Ay — AiV^^^ + \2V'''^\ where 

^(1) ^ ^.(H, + H^)y^^^^-i(H, + H^) ^ y(2) ^ + 26-^^^=+^.'). (53) 

Notice that F^^) (F^^^) does not depend on H^' {H^). Since V^^'> and F^^) 
uncorrelated quadratic averages separate as 

X-ly^^V,,) = \1{V^'V^') + Xlivl^V^P). (54) 
This leads to a natural separation of each of the contributions to purity 

l-{P{t)) = l- PW, it) + l- Pipl it) , (55) 

where Pspcc{t) denotes the average purity with particle i being a spectator, as given 
in section 14.11 In this way, the problem reduces to that of the spectator model. The 
respective expressions in section [4.11 mav be used. For instance, if we assume broken 
TRI, we obtain from ^ 



(P(t)) = l-4^A? fdr f\r'U]l+g^^lcosAy]a{r') + 0{X\N-'), (56) 
^^-^ Jo Jo ^ J 
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where Ci and C2 are the correlation functions of the corresponding environments 
defined in exact correspondence with psp . for and Hc> respectively. If in one or 
both of the qubits, the level splitting in the internal Hamiltonians is very large/small 
compared to the Heisenberg time in the corresponding environment (denoted by Tc 
and To' for and TJc', respectively) the degenerate and/or fast approximations may 
be used. As an example, if Ai ^ 1/tc and A2 <C 1/to' we find 

Fd(0 = 1 - (2 - 9e)i>4frSt) + A^/r„, W), (57) 
whereas if Ai ^ 1/tc and A2 ^ 1/tc' 



PAt) = 1 - [gi]lfrAt) + 2rc5S,*J - ^2 [gy/^Jr^, (t) + 2r.'9l%t\ ■ (58) 

It is interesting to note that if we have two separate but equivalent environments 
(i.e. both Heisenberg times are equal), we get exactly the same result as for a single 
environment. Also notice that the Hamiltonian of the entire system separates and 
thus the total entanglement of the two subsystems (Hi (g) He and H2 (E) He') becomes 
time independent. 

4-3. The joint environment configuration 

The last configuration we shall consider is the one of joint environment; see figureHJc). 
Its uncoupled Hamiltonian is 

Ho = Hi+H2 + He (59) 

whereas the coupling is given by 

AF = AiFe,l+A2K,2. (60) 

Notice the similarity with ([5T|l and ([5^ . However, as discussed in the introduction, 
they represent very different physical situations. The coupling in the interaction 
picture can again be split \V = AiF*^^^ + A^^-'V2, where 

1/(1) ^ Qi{Hi+H,)y^ ^g-i(//i+_ffe)^ y{2) ^ gi(//2+-ffo)y^ ^^-i{H2+H,)_ (^Qi^ 

Note the slight difference with ([55)1 . However still F'-"'^-' and F*^^-* are uncorrelated, 
enabling us to write again ([M]) . 

From now on, the calculation is formally the same as in the separate environment 
case. Hence we can inherit the result (|56|) directly, taking into account that since they 
come from the same environmental Hamiltonian, the two correlation functions are the 
same. In case any of the Hamiltonians fulfils the fast or degenerate limit conditions, 
the corresponding expressions to (j57p and (|58p can be written. As an example, if the 
first qubit has no internal Hamiltonian, and the second one has a big energy difference, 
the resulting expression for purity decay is 

{P{t)) = 1 - A? (2 - ge)frAt) - Xj [g<^J, /.„(*) + 2mgi^l^t] + OiX\ N'') , (62) 

where th is the Heisenberg time of the joint environment. Monte Carlo simulations 
showing the validity of the result were done with satisfactory results, comparable to 
those obtained in figure [Tj The parameter range checked was similar to that in the 
figure. 
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5. A relation between concurrence and purity 

In the previous section we iiave found tliat purity decay is very sensitive to the 
initial degree of entanglement between the two qubits. In other words, we related the 
behaviour of purity in time to the degree of entanglement (measured by concurrence) at 
t = 0. In this section, we study the relation between purity and entanglement (again 
measured by concurrence) as both are evolving in time. We focus on initial states 
which are Bell states such that purity and concurrence are both maximal and equal to 
one at t = 0. Since concurrence is defined in terms of the eigenvalues of a hermitian 
4 X 4-matrix, an analytical treatment even in linear response approximation is much 
more involved than in the case of purity. Instead, we work with a phenomenological 
relation between purity decay and concurrence discovered in [27j and further studied 
in [13]. It proves to be valid in a wide parameter range, and it allows to obtain an 
analytic prediction for concurrence decay. 

We study the relation between concurrence and purity on the CP-plane, where 
we plot concurrence against purity with time as a parameter. Since the initial state 
is a Bell state, the starting point is always at {C,P)t=o = (1,1). In this plane not 
all points are allowed for physical states; two constrains appear. The first one follows 
from the ranges of C G [0, 1] and P e [1/4, 1]. The second restriction follows from 
the "monogamy" of entanglement: if a qubit is very entangled with another qubit, 
it cannot be very entangled with the rest of the universe (i.e. the environment), 
implying some degree of purity. Conversely if the qubit is very entangled with the 
universe, it cannot be entangled with the other qubit. This statement can be made 
quantitative, which leads to the concept of maximally entangled mixed states. They 
define the upper bound of the set of admissible states (the corresponding area on the 
CP-plane is shown in figure [8l^ a) in grey). Another region of interest corresponds 
to those states, which form the image of a Bell state under the set of bi-local unital 
operations [28] (red area in the same figure). Finally, we have the Werner states 
pw = -I- (1 — a) I Bell) (Bell I, < a < 1, which define a smooth curve on the 
CP-plane (black solid line). The analytic form of this curve is [T5] 



and will be referred to as the Werner curve. Note that states belonging to the Werner 
curve are not necessarily Werner states. 

The spectator and separate environment configuration are strictly bi-local, 
whereas the joint environment configuration becomes bi-local in the large iVo limit. 
None of the three types of environmental couplings are strictly unital. In this limit 
we find from our numerical simulations that all CP-curves fall into the region defined 
by the bi-local unital operations. 

For figure El^b) , we perform numerical simulations of the spectator Hamiltonian 
assuming broken time reversal symmetry (GUE case). We compute the average purity 
(P(t)) as well as the average concurrence {C{t)) as a function of time. This average 
is over 15 different Hamiltonians and 15 different initial states for each Hamiltonian. 
We fix the level splitting in the coupled qubit to A = 1 and consider two different 
values A = 0.02 and 0.14 for the coupling to the environment, figure [8]^b) shows the 
resulting CP-curves for different dimensions of He- Observe that for both values of A, 
the curves converge to a certain limiting curve inside the red area (defined by the set 
of bi-local unital channels) as dim(?ic) tends to infinity. While for A ~ 0.02, this curve 




(63) 
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Figure 8. In (a) we show the area of concurrence-purity combinations which are 
allowed for physical states (grey area plus the set {{0, P),P G [1/4, 1/3]}). The 
image of a Bell state under the set of bi-local unital operations defines the red 
area. The Werner curve II63I I is shown as a thick black solid line. In (b) we show 
curves ((C(4)), {P{t))) as obtained from numerical simulations of the spectator 
Hamiltonian. GUE matrices are used during all numerical experiments in this 
section. We choose A = 1 and vary dim('Hc) for two different coupling strengths 
A = 0.02 (thin lines) and A = 0.14 (thick lines). The resulting curves are plotted 
with different colours, according to dim('He) as indicated in the figure legend. 



is at a finite distance of the Werner curve, for A = 0.14 it practically coincides with 
Cmv{P). To check this statement in more detail, consider the numerical CP-curve 
C'num(-P) obtained from our simulations and define its distance E to the Werner curve 
as 

dP|Cnum(P)-Cw(P)| . (64) 

The behaviour of E is shown in figure [Hlja). For A = 0.14 (black dots), the error 
goes to zero, which means that the corresponding CP-curves indeed converge to the 
Werner curve. In fact from a comparison with the black solid line we may conclude 
that the deviation E is inversely proportional to the dimension of Tie- By contrast, for 
A = 0.02 (red dots), the error tends to an approximately constant value, in line with 
the assertion that the numerical results converge to a different curve. In figurelSUb) we 
plot the error as a function of A, for different dimensions of Tie as indicated in the 
figure legend. The results suggest that the convergence to the Werner curve occurs as 
long as A > Ac which is of the order of 0.1 for A = 1. In fact for all 0.001 < A < 100 
studied, we found a critical value Ac such that the CP-curves converge to the Werner 
curve as long as A > Ac. For A = all studied couplings led to convergence to the 
Werner curve. 

In the presence of TRI the CP-curves again converge to the Werner curve, for 
values of the coupling greater than a critical value Ac- As in the GUE case, this critical 
value depends on A. However, in contrast to the GUE case, we find that Ac remains 
finite for A = 0. In the other configurations considered (the joint and the separate 
environment), the behaviour is similar. In those cases it is the largest (of the two) 
coupling strengths which determines whether the CP-curves converge to the Werner 
curve, or not. 

A partial explanation for this particular behaviour in the CP plane is provided 
in [55j. The authors find that for unital channels, the possible points in the CP plane 
that can be accessed are limited to a small 2-dimensional region, see figure ^a,) . 

We test the unitality of the channels induced by the non-unitary evolution in the 
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Figure 9. In (a) we show E defined in 1641 1 which measures the distance between 
the numerical CP-curves and the Werner curve as a function of dim A^c. The same 
data is used as in figure |8] The red dots refer to the case A = 0.02 where the error 
E apparently remains finite. The black dots refer to the case A = 0.14, where 
E seems to tend to zero. The black solid line, which is proportional to N^^ , is 
meant as a guide to the eye. In (b) we plot i? as a function of the coupling A, for 
various values of Ne- We may say that the line A 0.1 separates two regimes. 
For A > 0.1 we observe the convergence to the Werner curve, whereas for A < 0.1 
the limiting curve is a different one. 
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Figure 10. We evaluate the unitality condition for the spectator configuration 
(i.e. on one of the qubits). We vary the size of the environment and test different 
times, that would lead to different values of concurrence as shown in the inset, 
for the spectator case with the corresponding coupling. A line with slope —1/2 is 
also included for comparison. 



following way. Since we want to test = 1, we prepare an initial pure condition 
that leads to a completely mixed state in the qubit (let us focus on the spectator case) , 
i.e. , 



V2 



(65) 



with {ipe^^ ) = Sij. We let the state evolve with a particular member of the 
ensemble of Hamiltonians defined in Afterwards we evaluate the Euclidean 
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distance of the resulting mixed state in the qubit from the fully mixed state, within 
the Bloch sphere. The average distance is plotted as a function of the size of the 
environment in figure 1101 We conclude that the unitality condition is approached 
algebraically fast as the size of the environment increases. This by no means explains 
the generic accumulation to the Werner curve, but certainly restricts the possibilities; 
furthermore, for high purities/concurrences the area converges to the curve C = P, 
hence providing a satisfactory explanation in this regime. 
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Figure 11. We show numerical simulations for the average concurrence as a 
function of time in the joint environment configuration (GUE case). In (a) we 
consider small couplings Ai = A2 = 0.01 which lead to the Gaussian regime 
for purity decay. The red symbols show the result without internal dynamics 
(Ai = A2 = 0), whereas the black symbols are obtained for fast internal Rabi 
oscillations (Ai = A2 = 10). The corresponding theoretical expectation based on 
II67I I and II56I I is plotted as a red solid line (the case Ai = A2 = 0), and a black 
solid line (Ai = A2 = 10), respectively. In (b) we consider stronger couplings, 
Ai = A2 = 0.1, such that purity decay becomes essentially exponential (Fermi 
golden rule regime), while the level splitting have been set to Ai = A2 = 0.1 
(black symbols). The theoretical expectation (red solid line) is based on the same 
expressions as in (a). The insets in (a) and (b) display the corresponding evolution 
in the CP-plane with the same symbols as used in the main graph. In addition, 
the physically allowed region (grey area) and the Werner curve (black solid line) 
are shown. In all cases Nc = 1024 



Sufficiently close to P = 1, the above arguments imply a one to one 
correspondence between purity and concurrence, which simply reads C ~ P. This 
allows to write an approximate expression for the behaviour of concurrence as a 
function of time 

CLR(t) = PLR(t), (66) 

using the appropriate linear response result for the purity decay. The corresponding 
expressions have been discussed in detail in the previous section. ((66)) has similar 
limits of validity as the linear response result for the purity. Thus, in a sense, we may 
call it a linear response expression for concurrence decay. 

In those cases, where the coupling is beyond the critical coupling Ac and where 
the exponentiated linear response expression (j25p holds for the average purity, we can 
write down a phenomenological formula for concurrence decay, which is valid over the 
whole range of the decay 

CELR(t) - Cw(PELR(t)) . (67) 

In figure [IT] we show random matrix simulations for concurrence decay in the joint 
environment configuration. We consider small couplings Ai = A2 = 0.01 which lead to 
the Gaussian regime for purity decay, as well as strong couplings Ai = A2 = 0.1 which 
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lead to the Fermi golden rule regime. We find good agreement with the prediction 
(|67p. except for the Gaussian regime when we switch-on a fast internal dynamics in 
both qubits (Ai = A2 = 10). These level splitting are so large, that the coupling 
strengths are much smaller than Ac, which leads to the deviations observed. 

Note that, in the separate environment case, the entanglement of the two 
subsystems defined on the spaces Hi (E) Tic and 'H2 <8) Tic' is constant in time. For the 
initial conditions we use, the entanglement stems entirely from the two qubits. The 
concurrence of these two qubits decays because the constant entanglement spreads 
over all constituents of the two large subspaces. 

6. Conclusions and outlook 

We have derived general linear response formulae for the decay of purity for one and 
two qubits. The initial states are grouped in relevant families and the dependence on 
these families is analysed in detail. In particular we see that time reversal invariance 
will generally increase the number of families we have to consider, though entanglement 
will undo some of this in the two-qubit case. We do not discuss the well-known 
quadratic decay at very short (Zeno) times, though it obviously occurs in our model 
[6]. We find the linear decay resulting from the usual master equation treatment at 
short times, but we see the crossover to quadratic decay as we approach the Heisenberg 
time. This crossover is important if decoherence is dominated by a few degrees of 
freedom coupled to the central system, leading to a finite Heisenberg time in the 
relevant environment. We show that the quadratic term is absent if we choose an 
eigenstate of the intrinsic Hamiltonian of the central system. 

In general we can say that entanglement between the qubits accelerates 
decoherence mainly by affecting the prefactors in the linear response expression. 
Internal dynamics introduces an additional term but also affects the prefactors. It 
tends to stabilize the system though the interplay with entanglement does not make 
this statement strict for small changes. 

The analysis of concurrence decay and its relation to purity decay largely confirms 
previous findings |13ll27j in this more general setting. In particular the time-evolution 
follows the purity - concurrence relation of Werner states. We use this relation to 
give an analytic, albeit heuristic, formula for concurrence decay. Actually the relation 
shows small deviations for very weak couplings, and small concurrences. We obtain 
some understanding of these findings by showing that the process becomes unital in 
the limit of a large environment. 

We thus have given a rather complete description of the time-evolution of 
decoherence and entanglement of two non-interacting qubits. From a technical 
point of view we find a remarkable sequence of argumentation. In linear response 
approximation, all two-qubit configurations can be reduced to the spectator case. 
This case is equivalent to a single qubit, initially in a mixed state yet not entangled 
with the environment. The evolution of the mixed central system is obtained in the 
appendix in the linear response approximation. This result in turn allows to derive the 
general formulae for all cases mentioned above. While this allowed a simple unified 
representation of the two-qubit problem in the framework of an RMT model, actually 
one readily suspects, that this result is more general, and this is indeed the case as we 
shall see in a forthcoming paper [55]. Another line, that could be easily followed, is 
to consider a random interaction between the two qubits and a joint coupling to the 
environment. While such a study is easy and has some interest, the real problem to 
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tackle is to introduce time dependent operations such as gates, while the decoherence 
modelled by RMT acts. Another open problem that can be tackled in this framework 
is a situation of a near and far environment the first coupled weakly and the second 
extremely weakly. The first should probably have a finite Heisenberg time on a time 
scale relevant to the system, while the second more likely would have very dense 
spectrum and thus an infinite Heisenberg time. Finally, recall that we used formally 
echo dynamics for purely technical reasons. Yet this could be used to great advantage 
if we simultaneously want to study internal errors in the implementation of a quantum 
algorithm on a physical device as well as decoherence effects. As RMT is extremely 
successful in describing fidelity decay it seems to be the optimal approach. 
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Appendix A. One-qubit purity decay for pure and mixed states 

Here, we consider the one-qubit decoherence in a more general context. One- 
qubit decoherence as described in section [3] deals with a separable initial state 
go — Pi '9 |'0c)('0c|, where pi — |-0i)('0i| is a pure state of a single qubit coupled to 
an environment in the pure state IV'o). We consider here arbitrary, finite dimensions 
of Til and Tic, and we allow pi to be mixed. The latter poses some problems on the 
physical interpretation of purity as a measure for decoherence. However, these can be 
avoided by taking the point of view of the environment, as it becomes entangled with 
the central system. That means we consider the purity of the state of the environment 
after tracing out the central-system's degrees of freedom; see (|36|) . If we consider the 
entanglement with a spectator being responsible for the mixedness of pi , we can use 
the following results to describe purity decay in the spectator model. 

In this sense, the result can be applied to the spectator model described in 
section 14.11 

We will derive explicit expressions for the average purity {P{t)) as a function of 
time, where the average is only over the random coupling. These expressions involve 
a few elemental spectral correlation functions, whose properties will be discussed in 
[Appendix B[ As detailed in section 12.41 we work with the echo operator M (t) in 
the linear response approximation Note that I(t) is hermitian, while, due to 

time-order inversion, J{t) is not. 

The purity of the evolving state Let gi and g2 be two density matrices, defined on 
the Hilbert space TC' — Hi ^ Tic- We define the purity form p[-] as a function of pairs 
of such quantum states, 

p[gi ® g2\ = tro[tri(£ii) tri(£<2)] • (A.l) 

Since any linear operator acting on TC ® TL' can be expanded in terms of separable 
operators of the form gi® g2^ linearity implicitly defines the purity form for arbitrary 
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linear operators. For arbitrary operators A, B acting both on Ti.' , the purity form has 
the following property: 

p[A (g) B] = p[ylt B^]* = p[B (»A]= ^t]* (A.2) 

The purity defined in ([7|) can be expressed in terms of the purity form as 

P(tre q) = p[g (g) g] . (A.3) 

Note that on the RHS of this equation, we first take the trace over the central system, 
i. e. we consider the purity of the state of the environment after tracing out the central 
degrees of freedom. With this little twist we can describe the one-qubit decoherence 
in section [3] and the spectator model in section 14.11 at the same time. 

Purity decay in the linear response approximation To average purity, we use (jlip to 
compute g^^ (t) (E) g^'^ (t) in linear response approximation 

(P(t)) = {p[g^'it) ® g^'m ^p[{g''it) ® g^' (t))] . (A.4) 

Keeping terms only up to second order in A: 

g" (t) (g) g^Ht) = go®gQ 

-iA [/ go's/ go - go l'' ® go + go® I go - go® go l] ^ (A. 5) 

-A^ [j go® go + go ® go + go® J go + go® go J] ^ 

+a2 [lgoP®go-lgo®lgo + lgo®gol'' + gol''®lgo- go l'^ ® gol'' + go® I gol''] ■ 

In the next step, we perform the ensemble average over the coupling. As a result, the 
terms linear in A vanish. For the remaining terms, we use the properties in (|A.2p and 
the fact that I{t) is hermitian, to obtain 

(P(<)) = P(0) - A^ {Aj - Ai) P(0) - p[go ® go] 

Aj = 4Rep[(J) go ® go] 

Ai = 2{p[{I go I) ® go] - Rep[{I go® I go)] + p[{I go ® go I)] ) • (A.6) 

In order to pull out the same double integral of Aj and Aj, we use the time ordering 
symbol T. It allows to write for the average purity as a function of time 

(P(<)) = P(0) - 2A2 /" dT f dr'ReAji (A.7) 

An = p[T{VV')go ® go] - p[{V'goV) ® go] + p[{Vgo ® V' go)] 

~p[{V'go ® goV)] , (A.8) 
where V and V' are short forms for the coupling matrices V{t) and V{t'), respectively. 
The arguments t and t' of the coupling matrices are interchanged in the second and 
fourth term of Aji- This does not change the value of the integral of course, but it 
facilitates to handle common terms in the following considerations. 

The general result If the coupling matrix is taken either from the Gaussian unitary 
(GUE) or orthogonal (GOE) ensemble, we find 

An = [ Ci(|r - r'l) - S,ir - r') ] [C,{\t ~ r'|) - S,{r - r') ] 

+XGO^[l -S',{-t-t') + S[[t + t')S',[t + t')-S[{-t-t')] . (A.9) 

These expressions are derived in the following two subsections, for the GUE case 
in [Appendix A.l[ for the GOE case in [Appendix A.2[ The correlation functions 
Cx, Sx, C'^ and S'^ with a; G {1, e} are defined and discussed in [Appendix B 
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The limit dim(7ie) ^ oo; If the dimension of the Hilbert space of the environment 
goes to infinity, the corresponding correlation functions simphfy, as discussed 
in [Appendix B| We are then left with 

An= [Ci(|T-r'|)-5i(T-T')] C{\t-t'\) + xgoe[1-S[{-t-t')] .(A.IO) 
Appendix A.l. Appendix A.l The GUE case 

The integrand ylji in (|A.8[) consists of four terms. These terms will be considered, 
one after the other. We will always first average the argument of the purity form. 
Only then we perform the partial traces over subsystem Tii, and therefore the final 
trace over the environment. The averaging is only over the coupling, and it is done 
by applying two simple rules, as described below. 

For the coupling Vi.o in the product eigenbasis of Hq, we use either random 
Gaussian unitary (GUE) or orthogonal (GOE) matrices. Their statistical properties 
are completely characterized by the following second moments (for notational ease we 
ignore the subscript i.o for a moment) 

{VijVki) = 5ii Sjk + Xgoe Sik Sji , (A. 11) 

where xgoe = 1 if V" is taken from the GOE, while xgoe = if is taken from the 
GUE. In the interaction picture we then find: 

{V,j Vm) = ( S^l S,k e-'(^^-^') + XGOE 5,1 g-'^^^-^') ) . (A.12) 

p^T {VV') Qo ® Qo\{t)- We first compute the average 

T {V V' go ®go)^ ^ \ij) T (V^(t),j- fc; V{T')ki,rnn) {mn\ go (g) go 

ij kl ran 

= J2 Nj) e-'(^--^-) l^-^'l go(E>go, (A.13) 

ij kl 

where we have used the averaging rule, (|A.12[) . as well as the fact that the time- 
ordering operator T requests to exchange r with r' whenever r < t'. The indices i, k 
and m denote basis states in Hi, while the indices j, I and n denote basis states in He- 
We can rewrite that expression in a more compact form by employing the diagonal 
matrices Cx{t), defined in jBTTJ: 

T{VV' go(E>go) = (Ci(|T-T'|)®Ce(|r-T'|))(pi®|Vc) (^d )<^(pi®|^c) (^e|)-(A.14) 
Now we apply the partial traces over subsystem Hi 

tvi{ {VV')go) tn £<o = Ci(|t-t'|) Cc(|t - t'|) |i^c) (V-cI ■ (A.15) 
The final trace over the environment yields 

p[T{VV')go ® goW) - Ci(|t - t'|) CQt ~ r'\) . (A.16) 

p[{V' QoV) 8) Qo]{t): We first compute the average 

{V' goV)<Sigo^ I*-?') {%M {kl\go\mn) V^n,pq^ {pq\ ® go 

ij kl ran pq 

= E 1*^') e"'(^'^"^-) {kl\go\kl) ® go , (A.17) 

ij kl 



Decoherence of two qubit systems: A random matrix description 26 

where we have used (|A.12p . We apply the partial traces over subsystera Hi 
tniiV'goV)) tn eo = (E.fce''(^-^'=)(^-^') {k\pi\k)) 

X E (V'elO (Mo) (V'el 

ji 

= Ci(r-T') ^|j)e-(^^-^') ai^e) (V^elO (jI^c) (Vol . (A.18) 

The final trace over the environment yields 

piiV'goV) ® go] = Ci(t - r') ^(^^10 e'^' (--^') (/|Ve) (V-dj) e"'^^' (jlV-c) 

= Ci(r-T') S-olr-r') , (A.19) 
where S'a;(T) is defined in (|B.3p . 



^[(^£•0 (8i V'Qo)]{t): We first compute the average 
(Vgo^V'go) = J2 1*^') {^^^^^^' 

Po \mn) Vmn.pq^ {pq\ go 

ij kl run pq 

= J2 lu) e-'(^--^-) (fc^l eo ® \kl) go ■ (A.20) 

ij kl 

Then we apply the partial traces over subsystem TCi 

{tri{Vgo) tn {V'go)) = ( E.^e^^^''-^^) (^^^'^ (fc|pi|z) (^|pl|fc)) 

X E Ij) e-'(^'-^^) aiV-o) (V-cIO (Me) (V'el . (A.21) 

The final trace over the environment yields 

p[{Vgo ® V'go)] = S'i(r - r') S'e(r - r') . (A.22) 

p[{V'Qo ® QoV)]{t): We first compute the average 

{V' go^^goV) ^ 1*-^') (^4fei (^^1 ^0 «> |"^"-) V„,n^pq'^ {pq\ 

ij kl ran pq 

- E 1*^') e-'(^--^-H^-^') (fc^l go eo 1^0 (ul (A.23) 

We apply the partial traces over subsystem Hi 

{ tn {Vgo) tn (eo^') ) = (E,.e-'(^--^'=) (fc|pi|z) (^|pl|fc)) 

X E 1^') e-'(^^-^') m.) mi) 01 . (A.24) 

The final trace over the environment yields 

p[{Vgo ® goV')] = Si{t- t') Ce(r - r') . (A.25) 

Appendix A. 2. The GOE case 

In the GOE case, the average over the coupling yields, besides the previously 
considered GUE-term an additional one; see (jA.12p . In the following we will redo 
the calculation of the previous subsection, but consider only that additional term. As 
a reminder, we add the subscript 2 to the brackets which denote the ensemble average. 
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p[T{VV')2Qo ® Qo]{t): We first compute the average 

T{V V')2 Qo<^ 90= ^ T (V(T)y,fc; V{T)ki,„in)2 {mn\ go <^ go 

ij klmn 

= ^Vj){'ii\Qo®Qo = Qo®Qo, (A.26) 

where we have used the averaging rule, (jA.12p . Here the time-ordering operator has 
no effect. Now we apply the partial traces over subsystem Tii 

tri ( {V V')2 go ) tri go - (V'd ■ (A.27) 

The final trace over the environment yields 

p[{VV')2go®goKt) = l. (A.28) 

p[{V' QoV)2 ® Qo]{t): We first compute the average 

{V' go V)2 «) go = 1*-^') {^i'jM {kl\go\mn) V^n^pq^ ^ {pq\ «> go 

ij kl mn pq 

= J2 lu) e-'(^--^-) (^+^') {kl\go\tj) {kl\ ® go , (A.29) 

ij kl 

where we have used (jA.12p . We apply the partial traces over subsystem Tii 

tri{{V'goV)2) tri go = ^ |j) e^^^^-^') (^+^') (Zj^^e) (V^eb) (^l^c) (V-d .(A.30) 

ji 

The final trace over the environment yields 

pliV'goV) ® go] - ^'^''^ (^1^-^) (V'elj) e"'^^ (^+^') (^e|j) 

= 5^(-T-r'), (A.31) 
where S'^{t) is defined in (jB.4p . 



^[(^£•0 ® ^'^o)2](i)- We first compute the average 

{Vgo^V'go)2= J2 N-?') {^^^M {kl\ go ^ \mn) Vmn,pq^^ {pq\ go 

ij kl mn pq 

= J2 e-'(^"-^-) (^+^') {kl\ go ® {kl\ go . (A.32) 

ij kl 

Then we apply the partial traces over subsystem TCi 

{ in (Vgo) tri {V'go)), = ( E.^.e^'^^''-^-^ '^^^'^ {k\pi\i) (fc|pi|z)) 

X E 1^') (^+^') ai^c) (Velj) aiV'e) (V-el • (A.33) 

ji 

The final trace over the environment yields 

p[{Vgo ® V'go)2] = S[{t + r') S^r + r') . (A.34) 
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^[(^'£•0 ® QoV)2]{t): We first compute the average 

{V' go (E> go V)2 1*-^') (^4" (^^1 ® 1"^*^) Vrr,n,pq'^^ (mI 

ij kl mn pq 

= E 1*^') e"'(^-"^-) ® Nj) (fc/| . (A.35) 

ij kl 

We apply the partial traces over subsystem Hi 

{ tn {Vgo) tn {goV'))^ = { (^+^') (fclpi N) {k\p,\^) ) 

X E e-'(^^-^') (^+^') (^IV-c) (V'cb') . (A.36) 

The final trace over the environment yields 

p[{Vgo®goV')2]^S[i-T-T'). (A.37) 

Appendix B. Some particular correlation functions relevant for the 
average purity decay 

Appendix B.l. Definitions 

Averaging over the perturbation (i.e. the coupling) leads to expressions which may 
involve the diagonal matrix 

C.{t) = \i) e-'^^--'^^^ ^ {t\ x e 1, e . (B.l) 

ik 

Here, x denotes one of the two subsystems considered in [Appendix A[ either the qubit 
or the environment. Evidently, the energies are the eigenvalues of the corresponding 
Hamiltonian. In the derivations in [Appendix A[ the expectation value of Cx (t) with 
respect to the initial state go = pi®Pc are of particular importance. These are denoted 

by 

C,(t) =tr(C,(r) p,) a;el,e. (B.2) 

In this work, pc = |V'o)(V'o| is always a pure state. 

We will also encounter another type of correlation function, which may be defined 
as follows 

^x(t) =tr, [e-'^^>,e'^^>,] . (B.3) 

If the initial state px is pure, this quantity becomes the return probability. At r = it 
gives the purity of px- In the case of GOE averages, we also encounter the correlation 
function 

S'^{T)=iVx [e-'"-^Pxe'"-^pl,] . (B.4) 
Appendix B.2. Properties 

The limit of infinite dimension (of Tie): In the main part of this paper, we focus on 
the limit, where the dimension of the environment (s) becomes infinite. In that case 
it makes sense to perform an additional spectral average over He and/or Hgi. In the 
case of Cc{t) this yields 

(Ce(r)) = C(r) le C(r) = (Ce(T)) = l + 5(T/TH)-6f (t/th), (B.5) 
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where 62 (i) is the two-point spectral form factor with time measured in units of the 
Heisenberg time th for the corresponding spectral ensemble. In the limit of large 
dimension ~ dim(7ic) we find for the other correlation functions: 

5e(r) = tr (e-'^^'nV'c)(V'c|e'^"1Vc)(V'c|) (B.6) 

S',{r) = tr (e-'^-|V'e)(V'c|e'^-|Vc*)(^c1) ■ (B.7) 

For random pure states, as considered in the present work, both correlation functions 
are at most of order one at t = 0. For t > they drop very quickly and soon become 
of order N~'^. This happens on the same time scale, where Cc(t) drops from values 
of the order to values of the order one. In that sense we consider these correlation 
functions to contribute only 0{N^^) corrections to the result given in (|A.10p . 



The different correlation functions for a single qubit A single qubit is a two level 
system. The most general pure initial state is given by \a) ~ |l)ai + |2)a2, with 
|aiP + |a2p = 1- The most general mixed state is given by pi = Ai |a) (a| + A2 |6) {b\, 
with Ai,A2 > 0, real, Ai + A2 = 1, and |a), |6) arbitrary pure states with {a\b) = 
0. We will now investigate the behaviour of the different correlation functions 
ReCi(r), S'i(r), and S[{t) as it depends on the initial state and the Hamiltonian 

Fi = |i)ii;i(i| + |2)£;2(2|. 

(a) Assume pi = \a) {a\. Then 

CiiT)=J2\az\^e-''^^''~^'^" = l + cosAT + i(|a2|'- laiH sin At , (B.8) 

ik 

holds, so 

ReCi(T) = 1 + cosAr A ^ E2 - Ei . (B.9) 

(b) For pi = Ai \a) (a\ + A2 \b) (b\ we still have 

ReCi(T) = 1 + cosAr . (B.IO) 

(c) Assume pi = \a) {a\. Then we find for 5i(t) = s{\a);T) 

Siir) = J2 e-'(^'-^'='^ = |ai|4 + |a2|' + 2 \ai\^\a2\' cos Ar . (B.ll) 

ik 

This expression only depends on the absolute values squared of the coefficients 
ai and 02. Therefore we may parametrize them without loss of generality as 
ai = cos(/i and 02 = sini/). We then find 

•SiM =30 + (1-50) cos At ^ \ai\^ + \a2\^ = 1 - ^ sin^ 2(j) . (B.12) 

(d) For a general mixed state pi = Ai \a) (a| + A2 \b) {b\ we find 

5i(t) =A2s(|a);T)+A^s(|6);T)+2AiA2 ^a,<.6fe6* cos{E, - Ek)T 

ik 

= A2s(|a);T) +A^s(|6>;t) +4AiA2|ai|>i|2(l-cosAT) , (B.13) 

where we have used that {a\b) — aibl + 02^2 — 0- Note that the coefficients 
ai,a2,t'i,62 may be arranged into a square unitary matrix, and must therefore 
be of the following general form 

^ai bi\^,,( '''r.^, 'l""t) ■ (B.14) 
, ^2 02/ \ — e ^ sm e ^ cos r ' 
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This shows that s(|a); r) = s{\b); r) = + (1 — 5^) cos Ar, and that 

Si{t) = {Xl + Xl)[9^ + {l-g^) cosAr] + 2 A1A2 (1 - 30) ( 1 - cos At ) .(B.15) 

Since Ai + A2 = 1 it is convenient to set Ai — cos^ 9 and A2 — sin^ 9 such that we 
may write Aj + Aj = ge and obtain 

Si{T) = l^gg^g^ + 2ggg^ + {2gg-l){l-g^) cosAr. (B.16) 

(e) Assume again that pi — \a) {a\. Then we find for S[{t) — s'(|a);r) 

S[{t) = |ai|4 + |a2|4 + ajialf e'^^ + aHa^f e''^^ 
==50 + 2 |ai|2|a2p cos(At + 2t]) 

= 50 + (1 - 50) cos(At + 277) 77 = arg(ai) - arg(a2) . (B.17) 

However, as discussed in section 13.21 a natural symmetry around the y axis (in 
the Bloch sphere picture) appears for A = 0. In that case, using 7 as defined in 
that section, one can prove that 

1 - '^iM|a=o = (1 - - ^^0^(277)] = sin^ 7 (B.18) 
using elementary geometric and trigonometric considerations. 

(f) For a general mixed state pi ~ Xi \a) {a\ + A2 \b) {b\ we find 

S[{t) = a? s'{\a);T) + Xl s'(|6); r) + 2 A1A2 Re ^ a,46^6^ e^'^^'"^") ^ .(B.19) 

ik 

It follows fi'om (|B.14[) that arg(6i) — arg(&2) = arg(ai) — arg(a2) — tt, such that 
the equality s'(|a); r) = s'(|5); r) holds, just as in the case of Si{t). Therefore we 
may write 

S[{r) = ge[g^ + (1 - 50) cos(At + 277) ] + 2 A1A2 (1 - 50) Re( 1 + c'(^-+2"-) ) 
= 1-59-50 + 25950 + (259 - 1)(1 - 50) cos(Ar + 277) . (B.20) 

Note that the only difference to S{t) in case (d) is the additional phase 277 in the 
argument of the cosine function. Using the same angle 7 defined with (I45p . in 
complete analogy with (P7)l . and using (jB.lSp we obtain 

1 - S[ (t) I a=o = 1 - 5e + (25fl - 1) sin^ 7. (B.21) 
Appendix C. The exponentiation 

The linear response formulae derived throughout the article are expected to work for 
high purities or, equivalently, when the first two terms in the series (jlip approximate 
the echo operator well; it is of obvious interest to extend them to cover a larger range. 
The extension of fidelity linear response formulae has, in some cases, been done with 
some effort using super-symmetry techniques. This has been possible partly due to 
the simple structure of fidelity, but trying to use this approach for a more complicated 
object such as purity seems to be out of reach for the time being. Exponentiating 
the formulae obtained from the linear response formalism has proven to be in good 
agreement with the exact super-symmetric and/or numerical results for fidelity, if the 
perturbation is not too big. The exponentiation of the linear response formulae for 
purity can be compared with Monte Carlo simulations in order to prove its validity. 
We wish to explain the details required to implement this procedure in this appendix. 
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Given a linear response formula Plr(^) (for which P(0) = 1), and an expected 
asymptotic value for infinite time P^o, the exponentiation reads 



fELR(i) = Poo + (1 - Poo) exp 



1 



(C.l) 



This particular form guaranties that Pelr, equals Plr for short times, and that 
limj^oo PelrC^) = Poo The particular value of Poc will depend on the physical 
situation; in our case it will depend on the configuration and on the initial conditions. 
Let us write the initial condition as 

\^i2{0)) = cos0|OiO2) + sin^liiiz) (C.2) 

for some rotated qubits |0i), |li). In each of the qubits interacting with the 
environment wc assumed that for long enough times, the totally depolarizing channel 
f d will act (recall that the totally depolarizing channel is defined as £d [p] = 1/ tr 1 
for any density matrix p). Hence, for the spectator configuration, the final value of 
purity is 

Poo = P(^d ® l[|Vl2(^))(Vl2(^)|]) = f (C.3) 

whereas for both the joint and the separate environment configuration we use the 
value 

Poo = P (fd ® fd[|^12(^))(Vl2(^)|]) = l- (C.4) 

Good agreement is found with Monte Garlo simulations for moderate and strong 
couplings. 
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